Integrative Proteomics and Metabolomics Analysis Reveals the Role of Small Signaling Peptide Rapid Alkalinization Factor 34 (RALF34) in Cucumber Roots

The main role of RALF small signaling peptides was reported to be the alkalization control of the apoplast for improvement of nutrient absorption; however, the exact function of individual RALF peptides such as RALF34 remains unknown. The Arabidopsis RALF34 (AtRALF34) peptide was proposed to be part of the gene regulatory network of lateral root initiation. Cucumber is an excellent model for studying a special form of lateral root initiation taking place in the meristem of the parental root. We attempted to elucidate the role of the regulatory pathway in which RALF34 is a participant using cucumber transgenic hairy roots overexpressing CsRALF34 for comprehensive, integrated metabolomics and proteomics studies, focusing on the analysis of stress response markers. CsRALF34 overexpression resulted in the inhibition of root growth and regulation of cell proliferation, specifically in blocking the G2/M transition in cucumber roots. Based on these results, we propose that CsRALF34 is not part of the gene regulatory networks involved in the early steps of lateral root initiation. Instead, we suggest that CsRALF34 modulates ROS homeostasis and triggers the controlled production of hydroxyl radicals in root cells, possibly associated with intracellular signal transduction. Altogether, our results support the role of RALF peptides as ROS regulators.

The Arabidopsis RALF family consists of 37 members [15]. The Arabidopsis RALF34 (AtRALF34) peptide was suggested to be involved in the gene regulatory network controlling lateral root initiation [14]. AtRALF34, regulated by auxin and likely by ethylene, is expressed before the first asymmetric divisions of pericycle cells at the xylem pole [14],  (GUS control) and in roots with an overexpression construct (RALF34-OE). RT-qPCR analysis was performed using RNA isolated from the control group (mix of roots) and from individual transgenic roots. Statistical analysis using Wilcoxon's signed-rank test showed a significant increase (*, p < 0.001) of expression levels in the overexpression group compared with the control. The y-axis indicates the relative transcript level (2 -ΔΔCT method). (B) The Lateral Root Initiation Index in three outer cortical layers of control roots (GUS control) and of roots overexpressing RALF34 (RALF34-OE). Statistical analysis using Wilcoxon's signed-rank test showed no significant differences (p > 0.05) in the ILRI of individual roots in the overexpression group compared with control roots.

Treatment of Cucumber Roots with Synthetic CsRALF34 Affects Root Growth But Not Lateral Root Initiation
Wild-type cucumber roots were treated with 2 µM CsRALF34 peptide for 48 h ( Figure  2A). Before the start of the experiment, no statistically significant differences in root length were observed between cucumber seedlings in the control and in the peptide-treated group; the average initial root length was 6.03 cm for the control and 6.11 cm for the pep-  (GUS control) and in roots with an overexpression construct (RALF34-OE). RT-qPCR analysis was performed using RNA isolated from the control group (mix of roots) and from individual transgenic roots. Statistical analysis using Wilcoxon's signed-rank test showed a significant increase (*, p < 0.001) of expression levels in the overexpression group compared with the control. The y-axis indicates the relative transcript level (2 -∆∆CT method). (B) The Lateral Root Initiation Index in three outer cortical layers of control roots (GUS control) and of roots overexpressing RALF34 (RALF34-OE). Statistical analysis using Wilcoxon's signed-rank test showed no significant differences (p > 0.05) in the I LRI of individual roots in the overexpression group compared with control roots. Comparison of average root length between the control group (mock, n = 100) and treated roots (+CsRALF34, n = 100). Statistical analysis using Wilcoxon's signed-rank test showed no significant differences in root length before treatment (ns) and significant differences in root length (*, p < 0.001) after treatment. (B) Comparison of the Lateral Root Initiation Index in three outer cortical layers. Statistical analysis using Wilcoxon's signed-rank test showed no significant differences (p > 0.05) between the ILRI values of control roots (mock, n = 14) and treated roots (+CsRALF34; n = 10).

Overexpression of CsRALF34 Does Not Affect the Expression of CsGATA14, CsGATA24, and CsE2F/DP Genes
During the overexpression experiment (2.1), we also analyzed the expression levels of representatives of the CsE2F/DP family, i.e., CsE2Fa, etc., and of the cucumber orthologs of AtGATA23. The cucumber proteome contains six members of the E2F/DP protein family (CsE2F/DP) ( Figure 3). The expression of the six cucumber CsE2F/DP genes encoding these proteins were analyzed in control roots and in roots overexpressing CsRALF34. CsRALF34 overexpression did not affect the expression levels of any of these genes ( Figure 4). Expression of levels of the GATA genes, CsGATA14 and CsGATA24, previously proposed to represent AtGATA23 orthologs in cucumber [46], were examined as well. The expression levels of both genes did not change in response to CsRALF34 overexpression ( Figure 4). (A) Comparison of average root length between the control group (mock, n = 100) and treated roots (+CsRALF34, n = 100). Statistical analysis using Wilcoxon's signed-rank test showed no significant differences in root length before treatment (ns) and significant differences in root length (*, p < 0.001) after treatment. (B) Comparison of the Lateral Root Initiation Index in three outer cortical layers. Statistical analysis using Wilcoxon's signed-rank test showed no significant differences (p > 0.05) between the I LRI values of control roots (mock, n = 14) and treated roots (+CsRALF34; n = 10).

Overexpression of CsRALF34
Does Not Affect the Expression of CsGATA14, CsGATA24, and CsE2F/DP Genes During the overexpression experiment (2.1), we also analyzed the expression levels of representatives of the CsE2F/DP family, i.e., CsE2Fa, etc., and of the cucumber orthologs of AtGATA23. The cucumber proteome contains six members of the E2F/DP protein family (CsE2F/DP) ( Figure 3). The expression of the six cucumber CsE2F/DP genes encoding these proteins were analyzed in control roots and in roots overexpressing CsRALF34. CsRALF34 overexpression did not affect the expression levels of any of these genes ( Figure 4). Expression of levels of the GATA genes, CsGATA14 and CsGATA24, previously proposed to represent AtGATA23 orthologs in cucumber [46], were examined as well. The expression levels of both genes did not change in response to CsRALF34 overexpression ( Figure 4).    . Relative transcript levels of Cucumis sativus RALF34, GATA14, GATA24, and E2F/DP genes in control roots and in roots overexpressing CsRALF34. GUS control, transgenic roots expressing GUS; RALF34-OE, transgenic roots overexpressing CsRALF34. Statistical analysis using Wilcoxon's signedrank test showed significant differences (*, p < 0.05) between expression levels in the overexpression group compared with the control only for CsRALF34 (asterisk). The y-axis indicates relative transcript levels (2 −∆∆CT method).

Analysis of Biochemical Stress Response Markers in Cucumis Sativus Roots
The effects of the RALF34 peptide might be related to the plant stress response. To obtain a preliminary overview of the dynamics of lipid peroxidation, membrane integrity, and the status of antioxidant defense in control roots vs. roots overexpressing CsRALF34, several markers of the cellular antioxidant defense status, i.e., the tissue levels of hydrogen peroxide, lipid hydroperoxides, thiobarbituric acid-reactive substances (expressed as malondialdehyde equivalents, TBARS), and ascorbic acid/dehydroascorbate balance, were analyzed ( Figure 5). The results clearly indicated that hydrogen peroxide levels were significantly lower in roots overexpressing CsRALF34 in comparison to control roots ( Figure 5A). In contrast, the contents of thiobarbituric acid-reactive substances were increased almost two-fold in roots overexpressing CsRALF34 compared with control roots ( Figure 5B). On the other hand, the levels of lipid hydroperoxides and ascorbate were not affected by the overexpression of CsRALF34 ( Figure S1-2A,B). Interestingly, the contents of ascorbic acid were below the quantification limit in all samples, i.e., the ascorbate pool was represented only by the oxidized form of ascorbate, dehydroascorbate (data not shown). . Relative transcript levels of Cucumis sativus RALF34, GATA14, GATA24, and E2F/DP gen in control roots and in roots overexpressing CsRALF34. GUS control, transgenic roots expressi GUS; RALF34-OE, transgenic roots overexpressing CsRALF34. Statistical analysis using Wilcoxo signed-rank test showed significant differences (*, p < 0.05) between expression levels in the over pression group compared with the control only for CsRALF34 (asterisk). The y-axis indicates relat transcript levels (2 −ΔΔCT method).

Analysis of Biochemical Stress Response Markers in Cucumis Sativus Roots
The effects of the RALF34 peptide might be related to the plant stress response. obtain a preliminary overview of the dynamics of lipid peroxidation, membrane integri and the status of antioxidant defense in control roots vs. roots overexpressing CsRALF3 several markers of the cellular antioxidant defense status, i.e., the tissue levels of hydrog peroxide, lipid hydroperoxides, thiobarbituric acid-reactive substances (expressed malondialdehyde equivalents, TBARS), and ascorbic acid/dehydroascorbate balan were analyzed ( Figure 5). The results clearly indicated that hydrogen peroxide levels we significantly lower in roots overexpressing CsRALF34 in comparison to control roots (F ure 5A). In contrast, the contents of thiobarbituric acid-reactive substances were increas almost two-fold in roots overexpressing CsRALF34 compared with control roots (Figu 5B). On the other hand, the levels of lipid hydroperoxides and ascorbate were not affect by the overexpression of CsRALF34 ( Figure S1-2A,B). Interestingly, the contents of asco bic acid were below the quantification limit in all samples, i.e., the ascorbate pool w represented only by the oxidized form of ascorbate, dehydroascorbate (data not shown Figure 5. Biochemical characterization of Cucumis sativus transgenic control roots (p35S::gusA, GU control) compared with roots overexpressing CsRALF34 (p35S::CsRALF34, RALF34-OE). (A) Hyd gen peroxide contents in control roots and roots overexpressing CsRALF34; (B) contents of thiob bituric acid-reactive substance (TBARS) in control roots and roots overexpressing CsRALF34 (det mined as malondialdehyde equivalents). Asterisks denote statistically significant differences tween groups of samples, t-test: p < 0.05. The raw data are summarized in Supplementary Inf mation S2.

Analysis of Primary Metabolites
The analysis of primary metabolites relied on two methods: gas chromatograph quadrupole mass spectrometry with EI ionization (GC-EI-Q-MS) and ion-pair reverse phase high-performance liquid chromatography coupled online with triple quadrupo tandem mass spectrometry (RP-IP-HPLC-QqQ-MS/MS). After processing, both datas were merged into one matrix prior to statistical evaluation.
The first approach, GC-EI-Q-MS, aimed to identify and relatively quantify therma stable primary metabolites and revealed 321 total compounds, annotated in the aqueo Figure 5. Biochemical characterization of Cucumis sativus transgenic control roots (p35S::gusA, GUScontrol) compared with roots overexpressing CsRALF34 (p35S::CsRALF34, RALF34-OE). (A) Hydrogen peroxide contents in control roots and roots overexpressing CsRALF34; (B) contents of thiobarbituric acid-reactive substance (TBARS) in control roots and roots overexpressing CsRALF34 (determined as malondialdehyde equivalents). Asterisks denote statistically significant differences between groups of samples, t-test: p < 0.05. The raw data are summarized in Supplementary Information S2.

Analysis of Primary Metabolites
The analysis of primary metabolites relied on two methods: gas chromatographyquadrupole mass spectrometry with EI ionization (GC-EI-Q-MS) and ion-pair reversedphase high-performance liquid chromatography coupled online with triple quadrupole tandem mass spectrometry (RP-IP-HPLC-QqQ-MS/MS). After processing, both datasets were merged into one matrix prior to statistical evaluation.
The first approach, GC-EI-Q-MS, aimed to identify and relatively quantify thermally stable primary metabolites and revealed 321 total compounds, annotated in the aqueous methanolic extracts, obtained from cucumber control roots (GUS control) or from roots overexpressing CsRALF34 (RALF34-OE), respectively (Table S1-4). In this group, 234 individual compounds were identified as primary metabolites by spectral similarity search against available libraries and/or co-elution with authentic standards. Some metabolites appeared as several isomers and/or methoxyamine-trimethylsilyl derivatives; therefore, the total number of metabolites identified was only 175. The identified metabolites represented 27 amino acids, 6 amines, 15 fatty acids and esters, 29 organic acids (lactic, glycolic, salicylic, benzoic acids, as well as di-and tricarboxylic intermediates of the TCA cycle), 53 sugars (monosaccharides with uronic, aldonic acids, and phosphorylated derivatives, di-and oligosaccharides), 16 phenolic compounds, 3 heterocyclic compounds, 3 co-factors, 1 sterol, and 22 representatives of other classes.
After merging the two result sets, a combined matrix with 391 entries was built and processed using the MetaboAnalyst 5.0 online software tool. Hierarchical clustering analysis with a heatmap representation of the normalized relative abundances corresponding to the individual metabolites revealed essential intra-group variability in both the control group and in CsRALF34 overexpressing roots ( Figure 6A). Thus, variability within the control group was higher than in the CsRALF34-overexpressing group ( Figure S1-3). The principal component analysis did not show a complete separation of the two groups in the corresponding score plots, while 75.6% of the total variance could be explained by the first principal component (PC1, Figure 6B). against available libraries and/or co-elution with authentic standards. Some metabolites appeared as several isomers and/or methoxyamine-trimethylsilyl derivatives; therefore, the total number of metabolites identified was only 175. The identified metabolites represented 27 amino acids, 6 amines, 15 fatty acids and esters, 29 organic acids (lactic, glycolic, salicylic, benzoic acids, as well as di-and tricarboxylic intermediates of the TCA cycle), 53 sugars (monosaccharides with uronic, aldonic acids, and phosphorylated derivatives, diand oligosaccharides), 16 phenolic compounds, 3 heterocyclic compounds, 3 co-factors, 1 sterol, and 22 representatives of other classes.
After merging the two result sets, a combined matrix with 391 entries was built and processed using the MetaboAnalyst 5.0 online software tool. Hierarchical clustering analysis with a heatmap representation of the normalized relative abundances corresponding to the individual metabolites revealed essential intra-group variability in both the control group and in CsRALF34 overexpressing roots ( Figure 6A). Thus, variability within the control group was higher than in the CsRALF34-overexpressing group ( Figure S1-3). The principal component analysis did not show a complete separation of the two groups in the corresponding score plots, while 75.6% of the total variance could be explained by the first principal component (PC1, Figure 6B). Despite the essential intra-group variability, the control group and the CsRALF34 overexpressors demonstrated clear differences in their metabolic profiles. A t-test with Benjamini-Hochberg false discovery rate (FDR) correction at p ≤ 0.05 revealed that 26 metabolites were significantly downregulated in transgenic cucumber hairy roots overexpressing the CsRALF34 gene ( Figure S1-4). A fold change (FC) analysis (with an FC cut-off of FC ≥ 1.5) revealed 33 upregulated and 15 downregulated metabolites, demonstrating abundance dynamics associated with CsRALF34 overexpression (Figure S1-5). Finally, the combination of these two criteria yielded ten metabolites, three of which were downregulated, while seven were upregulated upon the overexpression of CsRALF34 (Table 1, Figure 6C). D-mannose, cGMP, and gluconic and/or galacturonic acid appeared to be most responsive to CsRALF34 overexpression and demonstrated clear down-regulation in the RALF34-OE roots. The group of upregulated metabolites included nicotinic acid, adenosine, glucose-1-phosphate, arginine, orotidine 5 -monophosphate, ornithine, and 6-phosphogluconic acid.

Analysis of Semi-Polar Secondary Metabolites
The analysis of semi-polar secondary metabolites revealed 486 and 374 compounds detected in negative and positive ion mode, respectively, which were organized in two data matrices according to their ionization mode. The principal component analysis (PCA) of both datasets combined revealed no group separation in the corresponding score plots ( Figure S1-6A,B). As can be seen from the results of hierarchical clustering with heat map representation ( Figure S1-6C,D), this observation might be explained by high intra-group variability, which was comparable to the inter-group dispersion. However, as can be seen from the corresponding volcano plots ( Figure S1-6E,F), the majority of the analytes only exhibited minimal (within 1.5-fold) non-significant inter-group differences. Although several analytes demonstrated very strong FC responses to CsRALF34 overexpression, the intensity of the corresponding levels varied strongly within the groups. Thus, in contrast to the primary metabolome, the secondary metabolome was not significantly affected by CsRALF34 overexpression.

Protein Isolation and Tryptic Digestion
To ensure the efficient extraction of root proteins and the maximal possible coverage of the Cucumis sativus proteome, we selected the phenol-based protein extraction method. Determination of the protein concentrations in the obtained isolates revealed extraction yields in the range of 0.234-0.463 mg/g fresh weight (Table S1-3). The assay precision was determined by SDS-PAGE with 5 µg of protein per lane ( Figure S1-7). Gel-based crossvalidation by the densitometric assessment of Coomassie-stained gels revealed overall lane densities of 1.53 × 10 6 -1.82 × 10 6 arbitrary units (AU, RSD = 5.8%, Table S1-6). This RSD value was sufficient for sample amount normalization and did not require additional recalculation. The signal patterns observed in the electropherograms were similar between lanes ( Figure S1-7A,B). Tryptic digestion of the obtained protein isolates was considered to be complete, as no bands could be detected ( Figure S1-7C,D) [47].

Annotation of Cucumis sativus Proteins
In total, 10,378 peptides were identified by MS/MS spectra in the whole dataset. Among them, 8561 and 8780 peptides were identified in control roots and in roots overexpressing CsRALF34, respectively (Figure S1-8A, Supplementary information S3, Table S3-1). Among these, 6983 peptides (67.3%) occurred in both groups, while 1598 (15.4%) and 1797 (17.3%) peptides were unique for either control or CsRALF34-overexpressing roots, respectively. Based on these identifications, 2302 possible proteins could be annotated (2147 and 2120 in control roots and in CsRALF34 overexpressors, respectively, Figure S1-8B, Table S3-2), which represented 1946 non-redundant proteins (i.e., protein groups): 1821 vs. 1802 in control roots vs. CsRALF34 overexpressors ( Figure S1-8C, Table S3-3). Among these, 1965 proteins (85.4%) were common for both groups, 182 (7.9%) and 155 (6.7%) proteins were unique for controls or CsRALF34 overexpressors, respectively (Figure S1-8B, Table S3-2). However, all these identifications relied exclusively on MS/MS spectra and did not consider inter-group alignments, which can result in the identification of some proteins as group-specific, although the corresponding peptide signal was detectable in the other group as well, but not fragmented. Fortunately, this issue can be solved by label-free quantification.

Label-Free Quantification
In the first step, we characterized the observed differences by principal component analysis applied to the whole dataset ( Figure 7). The results indicated low intra-group dispersion and a high degree of explained inter-group variance, i.e., more than 93% of group-specific proteins provided reliable group separation by principal component 1 (PC1). On the PC2 axis, no inter-group separation could be observed.
lane densities of 1.53 × 10 -1.82 × 10 arbitrary units (AU, RSD = 5.8%, Table S1-6). This RSD value was sufficient for sample amount normalization and did not require additional recalculation. The signal patterns observed in the electropherograms were similar between lanes ( Figure S1-7A,B). Tryptic digestion of the obtained protein isolates was considered to be complete, as no bands could be detected ( Figure S1-7C,D) [47].

Annotation of Cucumis sativus Proteins
In total, 10,378 peptides were identified by MS/MS spectra in the whole dataset. Among them, 8561 and 8780 peptides were identified in control roots and in roots overexpressing CsRALF34, respectively (Figure S1-8А, Supplementary information S3, Table  S3-1). Among these, 6983 peptides (67.3%) occurred in both groups, while 1598 (15.4%) and 1797 (17.3%) peptides were unique for either control or CsRALF34-overexpressing roots, respectively. Based on these identifications, 2302 possible proteins could be annotated (2147 and 2120 in control roots and in CsRALF34 overexpressors, respectively, Figure  S1-8B, Table S3-2), which represented 1946 non-redundant proteins (i.e., protein groups): 1821 vs. 1802 in control roots vs. CsRALF34 overexpressors ( Figure S1-8C, Table S3-3). Among these, 1965 proteins (85.4%) were common for both groups, 182 (7.9%) and 155 (6.7%) proteins were unique for controls or CsRALF34 overexpressors, respectively (Figure S1-8B, Table S3-2). However, all these identifications relied exclusively on MS/MS spectra and did not consider inter-group alignments, which can result in the identification of some proteins as group-specific, although the corresponding peptide signal was detectable in the other group as well, but not fragmented. Fortunately, this issue can be solved by label-free quantification.

Label-Free Quantification
In the first step, we characterized the observed differences by principal component analysis applied to the whole dataset ( Figure 7). The results indicated low intra-group dispersion and a high degree of explained inter-group variance, i.e., more than 93% of group-specific proteins provided reliable group separation by principal component 1 (PC1). On the PC2 axis, no inter-group separation could be observed. Label-free quantification revealed 208 differently expressed proteins (t-test: p ≤ 0.05, FC ≥ 1.5, Supplementary information S4) in cucumber hairy roots overexpressing CsRALF34 among the 2302 annotated proteins. Among the differentially expressed proteins, 92 were upregulated, and 116 were downregulated in CsRALF34 overexpressors compared with control roots (Table S4-1 and Table 2, respectively). The most strongly upregulated proteins were thioredoxin (17.8-fold), glutathione peroxidase (17.2-fold), carnitine operon protein CaiE (16.6-fold), and wound/stress protein (16.5-fold) ( Table 2). The 70 kDa heat shock protein (21.1-fold), ubiquitin-like domain-containing CTD phosphatase (18.7-fold), and phloem lectin (18.2-fold) demonstrated the most pronounced downregulation ( Table 2). PCA of the set of differentially expressed proteins confirmed the inter-group separation observed with the whole root proteome based on principal component 1 (PC1, 99.1% of the explained variance can be seen in the score plot presented in Figure S1-9A). To assess the contribution of individual proteins to the differences associated with CsRALF34 overexpression, the corresponding loading plot, and biplot were analyzed ( Figure S1-9B,C). The top 25 loadings which most strongly contributed to the CsRALF34-associated alterations in the cucumber root proteome are listed in Table 2 (all loadings of differentially expressed proteins are listed in Table S1-7). Among them, thioredoxin, glutathione peroxidase, acidic endochitinase, tripeptidyl peptidase II, and protein transport protein SEC23 were upregulated (FC range 17.8-1.9-fold) and 70 kDa heat shock protein, ubiquitin-like domain-containing CTD phosphatase, phloem lectin, 10 kDa chaperonin, NADH dehydrogenase, and UDP-glycosyltransferase 1 were downregulated (FC range 21.1-1.6-fold).

Functional Annotation of CsRALF34-Regulated Proteins
Functional annotation of the prospective CsRALF34-dependently regulated proteome (i.e., 208 proteins identified as differentially expressed) relied on Mercator4 software. The analysis revealed 26 functional classes (bins, Figure 8, Supplementary Information S5) among the 29 available. The up-or downregulated proteins were represented by 25 and 22 bins, respectively (i.e., proteins involved in the DNA damage response, multi-process regulation, and solute transport were downregulated, whereas the polypeptides related to chromatin organization, photosynthesis, RNA biosynthesis, secondary metabolism, and solute transport were upregulated). The analysis revealed 26 functional classes (bins, Figure 8, Supplementary Information S5) among the 29 available. The up-or downregulated proteins were represented by 25 and 22 bins, respectively (i.e., proteins involved in the DNA damage response, multi-process regulation, and solute transport were downregulated, whereas the polypeptides related to chromatin organization, photosynthesis, RNA biosynthesis, secondary metabolism, and solute transport were upregulated). The group of upregulated proteins was dominated by those involved in protein biosynthesis (15 entries), with ribosomal proteins L19, L18a, L35, S24, and eukaryotic translation initiation factor 3 showing the highest abundance gain in CsRALF34 overexpressors (17.0-, 16.4-, 15.7-, 16.7-, and 16.3-fold, respectively). Proteins involved in vesicle tracking also strongly contributed to the upregulated group with eight entries in total (for example, vacuolar-sorting receptor 7, exocyst complex component, and general vesicular transport factor p115 with 14.9-, 14.9-, or 13.1-fold changes in abundance) were part of this group. For 12 proteins (13.3% of all upregulated species), no function could be assigned in terms of the available bins; these were indicated as "not assigned" (Figure 8).
The group of downregulated polypeptides was dominated by those involved in protein homeostasis (15 entries) with 70 kDa heat shock protein, 10 kDa chaperonin, subtilisin-like serine protease, and Clp protease showing a decreased abundance in CsRALF34 overexpressors (21.1-, 17.6-, 17.4-, and 16.2-fold, respectively). Enzymes involved in protein modification also strongly contributed to the downregulated group, represented by seven entries. These included glycylpeptide N-tetradecanoyltransferase and dolichyl-diphosphooligosaccharide-protein glycosyltransferase with 16-and 15.4-fold abundance changes, respectively. For 22 downregulated proteins (17.6%), the function could not be established, and they were indicated as "not assigned" (Figure 8).   The group of upregulated proteins was dominated by those involved in protein biosynthesis (15 entries), with ribosomal proteins L19, L18a, L35, S24, and eukaryotic translation initiation factor 3 showing the highest abundance gain in CsRALF34 overexpressors (17.0-, 16.4-, 15.7-, 16.7-, and 16.3-fold, respectively). Proteins involved in vesicle tracking also strongly contributed to the upregulated group with eight entries in total (for example, vacuolar-sorting receptor 7, exocyst complex component, and general vesicular transport factor p115 with 14.9-, 14.9-, or 13.1-fold changes in abundance) were part of this group. For 12 proteins (13.3% of all upregulated species), no function could be assigned in terms of the available bins; these were indicated as "not assigned" (Figure 8).
The group of downregulated polypeptides was dominated by those involved in protein homeostasis (15 entries) with 70 kDa heat shock protein, 10 kDa chaperonin, subtilisinlike serine protease, and Clp protease showing a decreased abundance in CsRALF34 overexpressors (21.1-, 17.6-, 17.4-, and 16.2-fold, respectively). Enzymes involved in protein modification also strongly contributed to the downregulated group, represented by seven entries. These included glycylpeptide N-tetradecanoyltransferase and dolichyldiphosphooligosaccharide-protein glycosyltransferase with 16-and 15.4-fold abundance changes, respectively. For 22 downregulated proteins (17.6%), the function could not be established, and they were indicated as "not assigned" (Figure 8).
The prediction of subcellular localization relied on WoLF PSORT with subsequent manual verification based on database and literature data ( Supplementary Information S6). The results showed that the cytosol, plastids, and nucleus most strongly contributed to the CsRALF34-affected proteome (33.3%, 18.1%, and 17.6%, respectively ( Figure 9A,B). On the other hand, the cell wall, peroxisome, oil body, and cytoskeleton proteins (1.9%, 1.9%, 1.0%, and 0.5%, respectively) were least affected by CsRALF34 overexpression in cucumber roots. Interestingly, the localization patterns of up-and downregulated proteins were rather similar. The prediction of subcellular localization relied on WoLF PSORT with subsequent manual verification based on database and literature data ( Supplementary Information  S6). The results showed that the cytosol, plastids, and nucleus most strongly contributed to the CsRALF34-affected proteome (33.3%, 18.1%, and 17.6%, respectively ( Figure 9A,B). On the other hand, the cell wall, peroxisome, oil body, and cytoskeleton proteins (1.9%, 1.9%, 1.0%, and 0.5%, respectively) were least affected by CsRALF34 overexpression in cucumber roots. Interestingly, the localization patterns of up-and downregulated proteins were rather similar.

The Effects of CsRALF34 on Root Metabolic and Signaling Pathways
The metabolic and signaling pathways affected by CsRALF34 overexpression were examined by mapping individual regulated proteins using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Based on these data, the activation of specific signaling pathways, changes in sugar and tricarboxylic acid metabolism, and enhancement of protein biosynthesis and transport could be observed in root cells of CsRALF34 overexpressors compared with control roots (Figure 8, Supplementary Information S5).

The Effects of CsRALF34 on Root Metabolic and Signaling Pathways
The metabolic and signaling pathways affected by CsRALF34 overexpression were examined by mapping individual regulated proteins using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Based on these data, the activation of specific signaling pathways, changes in sugar and tricarboxylic acid metabolism, and enhancement of protein biosynthesis and transport could be observed in root cells of CsRALF34 overexpressors compared with control roots (Figure 8, Supplementary Information S5).

CsRALF34 Is Not Involved in the Lateral Root Initiation
In this study, an integrated multiomics approach was used to elucidate the role of the small signal peptide CsRALF34 in the development of cucumber root systems. The study was carried out on transgenic hairy roots of cucumber overexpressing CsRALF34 compared with transgenic hairy roots expressing p35S::gusA. In the roots overexpressing CsRALF34, as confirmed by RT-qPCR, the lateral root initiation index (ILRI) remained unchanged compared with transgenic control roots, indicating that the number of initiated lateral root primordia was not affected (Figure 1).
In contrast, the number of lateral root primordia was reduced in Arabidopsis ralf34 loss-of-function mutants [14]. Another effect on root development by AtRALF34 was a change in the root growth rate. A stunted root growth rate was described for plants treated with the synthetic AtRALF34 peptide compared with the control group [15]. Consistent with this result, Atralf34 mutant roots were characterized by an increase in growth rate compared with the roots of wild-type plants [51]. Individual developing hairy roots always differ from each other in age and growth rate; therefore, the effects of CsRALF34 overexpression on the growth rate of hairy roots cannot be determined. Therefore, wildtype cucumber seedlings' roots were treated with synthetic CsRALF34 to determine its effect on root growth. This effect was similar to that described for Arabidopsis [15]: treatment led to a statistically significant decrease in the root length increment (Figure 2). At the same time, the number of initiated lateral root primordia calculated using ILRI was unchanged in cucumber roots treated with the synthetic CsRALF34 peptide (Figure 2), as well as in transgenic roots overexpressing CsRALF34 (Figure 1).

CsRALF34 Is Not Involved in the Lateral Root Initiation
In this study, an integrated multiomics approach was used to elucidate the role of the small signal peptide CsRALF34 in the development of cucumber root systems. The study was carried out on transgenic hairy roots of cucumber overexpressing CsRALF34 compared with transgenic hairy roots expressing p35S::gusA. In the roots overexpressing CsRALF34, as confirmed by RT-qPCR, the lateral root initiation index (I LRI ) remained unchanged compared with transgenic control roots, indicating that the number of initiated lateral root primordia was not affected (Figure 1).
In contrast, the number of lateral root primordia was reduced in Arabidopsis ralf34 lossof-function mutants [14]. Another effect on root development by AtRALF34 was a change in the root growth rate. A stunted root growth rate was described for plants treated with the synthetic AtRALF34 peptide compared with the control group [15]. Consistent with this result, Atralf34 mutant roots were characterized by an increase in growth rate compared with the roots of wild-type plants [51]. Individual developing hairy roots always differ from each other in age and growth rate; therefore, the effects of CsRALF34 overexpression on the growth rate of hairy roots cannot be determined. Therefore, wild-type cucumber seedlings' roots were treated with synthetic CsRALF34 to determine its effect on root growth. This effect was similar to that described for Arabidopsis [15]: treatment led to a statistically significant decrease in the root length increment (Figure 2). At the same time, the number of initiated lateral root primordia calculated using I LRI was unchanged in cucumber roots treated with the synthetic CsRALF34 peptide (Figure 2), as well as in transgenic roots overexpressing CsRALF34 (Figure 1).
It has also been proposed that AtRALF34 might indirectly regulate the expression of AtGATA23 [14], the key marker gene of the early steps of lateral root initiation [22]. We previously proposed CsGATA14 and CsGATA24 to represent the putative cucumber orthologs of AtGATA23 [46]. Here, we show that neither CsGATA14 nor CsGATA24 expression levels were changed in response to CsRALF34 overexpression ( Figure 4).
Next, we tested the hypothesis that the expression of CsRALF34 may be involved in gene regulatory networks participating in the early steps of lateral root initiation. It is known that expression of the transcription factor E2Fa in Arabidopsis, a member of the E2F/DP family [52], is regulated by auxin-dependent activation of the transcription factors LBD18/LBD33 [53]. This is a common mechanism of the auxin-dependent activation of the mitotic cycle occurring in the xylem pole of the pericycle before the first asymmetric division during lateral root initiation. Phylogenetic analysis of Arabidopsis and cucumber E2F/DP (CsE2F/DP) proteins was carried out, leading to the identification of cucumber orthologs of Arabidopsis E2F/DP genes ( Figure 3). Expression levels of all CsE2F/DP genes were not changed in response to CsRALF34 overexpression ( Figure 4).
Taken together, the data obtained using I LRI calculations for transgenic hairy roots overexpressing CsRALF34, as well as for CsRALF34-treated roots, and the evaluation of CsE2F/DP, CsGATA14, and CsGATA24 expression suggest that CsRALF34 does not represent a key player in lateral root initiation occurring in the root apical meristem of the cucumber parental root. However, although CsRALF34 overexpression does not affect lateral root initiation (i.e., it does not change the number of initiated lateral roots), an effect on the root growth rate cannot be excluded at this point.
The absence of any visible phenotype in roots overexpressing CsRALF34 prompted us to perform a comprehensive analysis of changes in the proteome and metabolome in response to CsRALF34 overexpression. Protein isolation followed by tryptic digestion was performed. The variability in protein yields (0.234-0.463 mg/g fresh weight, Table S1-6) was relatively low, in agreement with the low amounts of protein in roots compared with other parts of the plant [54,55]. Importantly, the precise estimation of protein extraction recoveries (RSD = 5.8%) were in line with our prior results [56]. Therefore, the observed lower variability in protein yields and precise estimation of the protein concentrations in the extracts enabled the reliable implementation of label-free quantification in the analysis of heterogeneous root proteomes [47].
PCA applied to the whole proteome ( Figure 7) exhibited low intra-group variance and a high degree of inter-group variance (more than 93%). These results indicated a high likelihood that the observed differences were predominantly associated with the effects of CsRALF34 overexpression.

Activation of Protein Biosynthesis
Our experiments clearly demonstrated that the overexpression of CsRALF34 resulted in a pronounced upregulation of proteins involved in protein biosynthesis and transport ( Figure 12). Upregulation of proteins constituting the small and large ribosomal subunits, RNA-binding family protein, nucleolar protein 4, ribosome export adaptor, protein transporter Sec23/24, and others were observed ( Figure 12, Table S4-1).

Inhibition of Root Growth and Regulation of Cell Proliferation
To date, in several studied plant species (e.g., Medicago truncatula and Solanum lycopersicum), the main role of the RALF34 peptide was reported to be the control of the alkalization of the environment for growth inhibition by suppressing the functions of H + -ATPase [11,16,57]. This should also prevent cell elongation. However, in other studies, this effect was discussed as a result of the molecular signaling cascade activated by RALF34. Pearce, Moura, Stratmann, and Ryan [11] showed that the expression of RALF34 was as-sociated with a decrease in the activity of the mitogen-activated protein kinase (MAPK), which was associated with the inhibition of lateral root growth. The data presented here revealed that the overexpression of CsRALF34 in cucumber roots resulted in a decrease in CDKA levels, which, in turn, could lead to a G2/mitosis transition block [48]. ethylene in plant tissues combined with a dwarf phenotype. Therefore, the autho sumed that the FER receptor kinase interacts with S-adenosylmethionine synthase downregulates ethylene biosynthesis in response to environmental stress, as well a external application of auxin and brassinosteroids [70]. Remarkably, Deslaurier Larsen [71] showed that mutation of fer led to a loss of sensitivity to brassinosteroid resulted in enhancement of the ethylene response. Finally, Bergonci et al. [13] showe RALF1 can compete with brassinosteroids for components shared by both signal duction pathways.
In summary, the available data are consistent with our results. It is likely that binding to its receptor, RALF34 suppresses the expression of aminocyclopropane ca ylate synthase, which then leads to a decrease in ethylene levels.

Plant Material and Bacterial Strains
Cucumber (Cucumis sativus L.) cv. Phoenix (Sortsemovosch, Saint Petersburg, R was used in this study. Escherichia coli strain XL-1 Blue was used for molecular cloning. Rhizobium rhizo (Agrobacterium rhizogenes) strain R1000 was used for the genetic transformation of p

Impact of CsRALF34 on ROS Signaling and Stress Adaptation
The exact mechanisms behind the contribution of RALF peptides to stress response and adaptation are still unknown. However, Stegmann et al. [60] showed that several AtRALFs are involved in the regulation of intracellular ROS levels. The authors reported that AtRALF23 mediated an increase in ROS production, whereas AtRALF17, in contrast, caused a decrease in ROS levels. Our results generally support the role of RALF peptides as ROS regulators, suggesting that CsRALF34 modulates ROS homeostasis. At least CsRALF34 caused alterations in the function of the respiratory chain combined with activation of the antioxidant defense systems.
We observed decreases in the relative abundances of NADH-dehydrogenase, complex I of the respiratory chain, and succinyl-CoA synthetase, the enzyme that converts succinyl-CoA to succinate (Supplementary Information S4, Table S4-2). Succinate deficiency is known to suppress the activity of succinate dehydrogenase, which is an enzyme of the TCA cycle and simultaneously represents complex II of the mitochondrial respiratory chain. The loss of its activity leads to a deficit of electrons in the respiratory chain. Moreover, as NADH-dehydrogenase reduces ubiquinone, which, in its reduced form, is the substrate of cytochrome c reductase, the activity of complex III could also be negatively affected by CsRALF34. Thus, given that complexes I and III are the main generators of the superoxide anion radical [61,62], this decrease in the electron flow through the entire respiratory chain might lead to a decrease in the rate of superoxide generation.
As shown here, the levels of H 2 O 2 were significantly lower in roots overexpressing CsRALF34 compared with control roots ( Figure 5A). This can be explained by the inhibition of electron transfer along the respiratory chain described above because superoxide anions are converted to H 2 O 2 by superoxide dismutase. On the other hand, the increased levels of thioredoxin 1, heat shock proteins Hsp90 and Hsp72, ubiquitin carboxyl-terminal hydrolase 1, and universal stress protein indicate that the cellular antioxidant defense was upregulated.
Equally importantly, glutathione peroxidase levels were also upregulated in response to CsRALF34 overexpression (Table S4-1). This enzyme is directly involved in antioxidant defense, specifically in the reduction of lipid hydroperoxides to the corresponding alcohols and in the reduction of hydrogen peroxide to water [63]. Thus, it can be concluded that CsRALF34 is critically involved in the suppression of H 2 O 2 production and the activation of multiple cellular antioxidant systems. It can be assumed that CsRALF34 mediated these effects via signaling cascades.
Due to their upregulation upon CsRALF34 overexpression, we could identify protein phosphatase 2C (PP2C) and calcium-dependent protein kinase 12 (CDPK12) as the components of the signaling pathways involved in the RALF34-dependent responses. PP2C is involved in abscisic acid (ABA) signal transduction. Binding of ABA to its receptors, the pyrabactin resistance/pyrabactin resistance-like/regulatory components of the abscisic acid receptor (PYR/PYL/RCAR), leads the receptors to bind to PP2C, which, in turn, leads to its inhibition and the activation of the SNF1-related protein kinase 2 (SnRK2). In the absence of ABA, SnRK2 proteins are inhibited by PP2C-dependent dephosphorylation (Hsu et al., 2021). Briefly, PP2C plays an important role in ABA signaling [64]. Importantly, CDPK12 is known to contribute to regulating the expression of antioxidant defense genes; therefore, its activity might affect plant adaptations to stress [50]. Thus, the RALF34induced upregulation of PP2C could affect ABA signaling and, thereby, plant adaptation to osmotic stress. In this context, it is possible that the upregulation of the cellular antioxidant pathways, the decrease in ROS production, and other metabolic shifts (Figures 10 and 11) could be mediated by increases in PP2C expression levels.
The involvement of ROS in this process was supported by the analysis of the tissue oxidative status in control and CsRALF34 overexpressing roots ( Figures 5 and S1-2). Consistent with the upregulation of the ROS-detoxifying enzymes, a clear decrease in H 2 O 2 content was observed in CsRALF34 overexpressing roots. On the other hand, levels of thiobarbituric acid-reactive substances were increased in CsRALF34 overexpressors compared with wild-type roots, indicating a complex effect of RALF34 on different aspects of antioxidant defense.
This conclusion is in agreement with the work of Song et al. [65], who showed that a mutation in FERONIA (FER), the RALF receptor gene in Arabidopsis, was associated with lower levels of ROS formation. Thus, it appears to be likely that RALF might directly affect ROS production. Moreover, because CDPK is activated by interaction with Ca 2+ ions [66], upregulation of the production of this kinase might result in increased sensitivity of the root to intracellular calcium levels. Upon treatment with nanomolar concentrations of natural or synthetic RALF peptides, the increasing levels of cytoplasmic Ca 2+ reached their maximum within 40 s [67]. This permits us to speculate that different RALF peptides can act in a concerted way; some RALF peptides might bind to the RALF receptor and initiate the release of Ca 2+ , whereas the others could increase the levels of sensitivity to ABA. Finally, it is important to mention that Ca 2+ levels might fluctuate strongly even over short periods of time. Interference of these dynamically altering calcium levels with RALF-mediated responses might underlie a plethora of RALF effects.

RALF-Related Dynamics of Cellular Metabolism
Based on the results of our study, CsRALF34 appeared to affect the energy metabolism of the root cell. As mentioned above, the overexpression of CsRALF34 was accompanied by the downregulation of NADH-dehydrogenase and succinyl-CoA synthetase. It appears likely that a compensatory mechanism of ATP generation was in effect, related to the upregulation of the TCA cycle, glycolysis, pentose phosphate pathway, and purine biosynthesis de novo ( Figure 11); the effect of CsRALF34 treatment on the root growth increment was not dramatic.
Altogether, the analysis of the primary and secondary metabolome showed that the primary metabolome was much more affected by CsRALF34 overexpression ( Table 1). The levels of D-mannose and cGMP were decreased, whereas the levels of glucose-1-phosphate, 6-phosphogluconic acid, ornithine, and adenosine were upregulated (Table 1). These results are consistent with the proteome analysis data with regard to enzymes involved in sugar metabolism.

RALF34 as a Modulator of Phytohormone Responses
Currently, the effect of RALF peptides on hormonal regulation is discussed in the literature [12]. In the most comprehensive way, this aspect was addressed in the studies of fer mutants, i.e., the plants defect at the gene of receptor-like kinase FER, which is one of the RALF receptors [28]. As was shown recently, this receptor is involved in the modulation of jasmonic acid, ethylene, ABA, and brassinosteroid dynamics in plants [28].
The results of this study (Section 2.10, Figure 10) clearly indicate an upregulation of PPC2, which can be expected to lead to a downregulation of SnRK2. This, in turn, would affect responses to ABA, ethylene, jasmonic acid, and to abiotic stress in general ( Figure 10). This is consistent with the observation that a fer loss-of-function mutant demonstrated hypersensitivity to both ABA and abiotic stresses associated with exposure to high salt concentrations and low temperatures [68]. The results of this study support the assumption that RALF34 is involved in the response to abiotic stress.
In addition to the effects related to the ABA signaling pathway, the overexpression of CsRALF34 was associated with a downregulation of aminocyclopropane carboxylate oxidase, which is responsible for the final step of ethylene biosynthesis [69]. This fact might suggest a decrease in the ethylene levels in the roots overexpressing CsRALF34. Notably, Mao et al. [70] showed that Arabidopsis fer mutants displayed higher levels of ethylene in plant tissues combined with a dwarf phenotype. Therefore, the authors assumed that the FER receptor kinase interacts with S-adenosylmethionine synthases and downregulates ethylene biosynthesis in response to environmental stress, as well as the external application of auxin and brassinosteroids [70]. Remarkably, Deslauriers and Larsen [71] showed that mutation of fer led to a loss of sensitivity to brassinosteroids and resulted in enhancement of the ethylene response. Finally, Bergonci et al. [13] showed that RALF1 can compete with brassinosteroids for components shared by both signal transduction pathways.
In summary, the available data are consistent with our results. It is likely that after binding to its receptor, RALF34 suppresses the expression of aminocyclopropane carboxylate synthase, which then leads to a decrease in ethylene levels.

Plant Material and Bacterial Strains
Cucumber (Cucumis sativus L.) cv. Phoenix (Sortsemovosch, Saint Petersburg, Russia) was used in this study.
Escherichia coli strain XL-1 Blue was used for molecular cloning. Rhizobium rhizogenes (Agrobacterium rhizogenes) strain R1000 was used for the genetic transformation of plants.

Reagents
The full list of the reagents is given in Supplementary Information S1 (Protocol S1-1).

Phylogeny and Bioinformatics
Sequences of six Arabidopsis E2F proteins [72] and two Arabidopsis DP proteins [73] were downloaded from the Arabidopsis Information Resource (TAIR, www.arabidopsis.org, accessed on 10 March 2023) [74] and used as a query to find amino acid sequences of C. sativus (cucumber, cv. Chinese Long v2) [75] in the Cucurbit Genomics Database v1 (CuGenDB v1, cucurbitgenomics.org, accessed on 10 March 2023) [76]. All alignments were performed using online Clustal Omega software (www.ebi.ac.uk/Tools/msa/clustalo/, accessed on 10 March 2023) [77] at default settings. The alignment file was transferred into MEGA7.0 software [78], followed by phylogenetic tree construction using the Maximum Likelihood method [79]. The Jones-Taylor-Thornton model [80] with evolutionary rate differences among sites (+G parameter) (Yang, 1994) was used for phylogeny reconstruction of the E2F/DP family in cucumber. Phylogeny was tested using the bootstrap method with 1000 replicates.

RT-qPCR Assays
Total RNA was extracted from frozen plant material using ExtractRNA reagent (Evrogen, Moscow, Russia). RNA quantity and integrity were measured using a Qubit 4.0 fluorimeter (Thermo Fisher Scientific, Waltham, MA, USA) using Qubit RNA BR Assay and IQ Assay Kits. Reverse transcription was performed as described previously [46]. For each qPCR assay, 1 µL of cDNA (80-120 ng) from a non-diluted sample (total volume 20 µL) was used.
The RT-qPCR analysis was performed using a Quant Studio 5 Real-Time PCR system (Thermo Fisher Scientific) in a total volume of 20 µL. qPCR and PCR conditions were described previously [46]. Primers used for qPCR (Table S1-1) were designed using Vector NTI Advance v 11.0 software (Thermo Fisher Scientific). Purified PCR primers were purchased from Evrogen (Moscow, Russia). Quantification cycles (Cq) were determined using Quant Studio Design and Analysis software v. 1.5.1 (Thermo Fisher Scientific). Relative transcript levels were calculated as described previously [46]. Elongation factor EF1α was chosen as a reference gene according to data on the stability of reference gene expression in cucumbers presented in the literature [81].
RT-qPCR analysis of CsRALF34 relative expression levels in the first replicate of the overexpression experiment was performed five times independently from the control group (GUS control, mixed root sample) and for each individual transgenic root (RALF34-OE, n = 9). The second replicate of the CsRALF34 overexpression experiment was conducted to analyze the transcript levels of CsRALF34 as well as the CsGATA14, CsGATA24, and CsE2F/DP genes in response to CsRALF34 overexpression. RT-qPCR analysis of the second replicate was performed four times independently.

Lateral Root Initiation Index (I LRI ) Estimation
The Lateral Root Initiation Index (I LRI ) [82] was estimated for two groups of roots: individual transgenic roots overexpressing CsRALF34 (n = 8) and GUS control roots (n = 8), as well as for wild-type roots treated with synthetic CsRALF34 peptide (n = 10) and for corresponding wild-type control roots (n = 14). The I LRI was calculated on longitudinal root sections as described previously [83]. The most proximal lateral root primordia were detected in the region located 2-5 mm from the root tip above the division zone of the cortex [84].

Treatments with Synthetic CsRALF34 Peptide
A synthetic CsRALF34 peptide (FWRRVHYYISYGALSANRIPCPPRSGRPYYTHN-CYKARGPVNPYTRGCSAITRCRR; >87% purity) was synthesized by ProteoGenix (Schiltigheim, France), and its structure was validated independently by mass spectrometry. Lyophilized peptide was diluted in a mixture of acetonitrile and double distilled water (1:3). Five-day-old wild-type cucumber seedlings with 5-7 cm long roots were incubated in aerated Hoagland's medium [85] supplemented with 2 µM CsRALF34 [15,86,87] or with the same volume of the solvent (control) for 48 h. Each experiment included 50 seedlings (both in the control and in the group with CsRALF34 treatment) and was repeated twice independently. For each seedling, the root length was measured before and after peptide treatment.

Molecular Cloning and Vector Design
Two genetic constructs were developed for overexpression assays using multisite Gateway technology (Gateway LR Clonase II plus, Thermo Fisher Scientific): for CsRALF34 overexpression using the 35SCaMV promoter and a control vector for gusA (β-glucuronidase A) overexpression under the control of the same promoter. The pKGW-RR-MGW binary vector, carrying the pAtUBQ10::DsRED1 [88] screening cassette in the backbone, was used as the destination vector (kindly provided by Erik Limpens, Wageningen University, Wageningen, The Netherlands). LR plus clonase reactions were prepared according to the manufacturer's instructions. The sequences of the cassettes of all constructs were verified by the PCR amplification of fragments and sequencing of the products. All primer sequences and their combinations are listed in Tables S1-2 and S1-3.

Plant Transformation and Handling for CsRALF34 Overexpression Assays
Agrobacterium rhizogenes (Rhizobium rhizogenes)-mediated transformation of cucumber seedlings was carried out as described previously [46,90] with minor modifications. The experiment was performed in two replicates. At least 25 composite plants with roots overexpressing CsRALF34 (p35S::CsRALF34) and GUS control plants (p35S::gusA) were obtained in each transformation. Transformants were cultured in vermiculite moistened with 4x Hoagland's medium. Transgenic roots 6-10 cm in length were harvested from both groups four times at 7-day intervals and flash-frozen in liquid nitrogen. The total number of roots was as follows: 220 for GUS control plants and 246 for CsRALF34 overexpressors.

Determination of Hydrogen Peroxide Contents
Analysis of H 2 O 2 contents in cucumber roots relied on the method described by Chantseva et al. [91], with modifications explained in Supplementary Information S1 (Protocol S1-2). The absorbance was measured at 575 nm (length of the optical path: 1 cm), and calculations were performed as described previously [91].

Determination of Lipid Peroxidation Product Contents
Lipid peroxidation products were quantified as malondialdehyde equivalents according to the protocol described by Soboleva et al. [92], with modifications presented in Supplementary Information S1 (Protocol S1-3).

Determination of Lipid Hydroperoxide Contents
Analysis of lipid hydroperoxide contents relied on the method of Frolov et al. [56], with modifications described in Supplementary Information S1 (Protocol S1-4).

Determination of Ascorbic Acid Contents
The determination of ascorbate contents and the respective statistical calculations were performed as described by Shumilina et al. [54] with modifications presented in Supplementary Information S1 (Protocol S1-5).

Protein Isolation and Determination
For protein extraction, approximately 250 mg of milled frozen root material was used, as described previously [95], with some modifications (Supplementary Information S1 (Protocol S1-6)). The precision of the assay was cross-verified by SDS-PAGE, according to Greifenhagen et al. [96], with some modifications. Briefly, after the gels were stained with 0.1% Coomassie Brilliant Blue G-250 for 12 h, the average densities across individual lanes (expressed in arbitrary units) were determined using a ChemiDoc XRS imaging system controlled by Quantity One 1-D analysis software (Bio-Rad Laboratories, Hercules, CA, USA). To calculate the relative standard deviations (RSDs), the densities of individual lines were normalized to the gel average value by ImageJ software [97].

Tryptic Digestion and Sample Pre-Cleaning
The day before the digestion procedure, Amicon Ultra 30K filter units (Sigma-Aldrich, Saint-Louis, MO, USA) were passivated in 500 µL 5% Tween-20 by shaking (350 rpm, 25 • C, 12 h). The next day, the filter units were washed twice with Milli-Q water while shaking (450 rpm, 30 min, 25 • C). Subsequently, 35 µg of protein was transferred to filter units, and 200 µL of urea solution (8 M urea, 50 mM Tris-HCl buffer, pH 7.5, UA) was added. Then the tubes with filters were centrifuged (14,000× g, 10 min, 4 • C), and the supernatants were discarded. All these steps with UA were repeated three times. Furthermore, 100 µL of reducing solution (100 mmol/L dithiothreitol in UA) was added, and incubated under continuous shaking (450 rpm, 25 • C, 1 h). Then, the mixtures were centrifuged (14,000× g, 10 min, 25 • C), and the supernatants were discarded. The same centrifugation and discarding steps were performed after the subsequent addition of 100 µL of alkylation solution (50 mmol/L iodoacetamide in UA); the incubation under continuous shaking (450 rpm, 25 • C, 1 h) was performed in the dark. The next step was washing filters with 200 µL of UA twice (centrifugation at 14,000× g, 10 min, 4 • C and discarding of supernatant after each step). Then, the same procedures were carried out with the digestion buffer (50 mmol/L ammonium bicarbonate, ABC) twice. The stock solution of trypsin (0.5 µg/µL) was freshly prepared, added to the solutions on the filters, and then incubated under continuous shaking (450 rpm, 37 • C, 4 h). The incubations were further continued after the addition of trypsin stock solution (enzyme-to protein-ratio w/w 1:50) under continuous shaking (450 rpm, 37 • C) overnight. The samples were centrifuged (14,000× g, 10 min, 25 • C), and the filtrates containing the mixtures of proteolytic peptides were supplemented with 40 µL of ABC buffer and centrifuged again (14,000× g, 10 min, 25 • C). This step was repeated twice; then, proteolytic hydrolysates were frozen at -20 • C. The completeness of the digestion was verified by SDS-PAGE [96].

Solid Phase Extraction
The proteolytic hydrolysates were pre-cleaned by reverse-phase solid-phase extraction (RP-SPE) using the elution scheme of Spiller et al. [98], with minor modifications specified in Supplementary Information S1 (Protocol S1-7).

Nano LC-MS/MS
The protein hydrolysates were loaded on an Acclaim PepMap 5 mm Trap Cartridge (Thermo Fisher Scientific) and separated on a Bruker FORTY separation column (C18 ReproSil AQ, 40 cm × 75 µm, 1.9 µm, 120 A; Bruker Daltonics, Bremen, Germany) using a nanoElute UHPLC chromatography system (Bruker Daltonics) coupled on-line to a TimsToF Pro quadrupole time-of-flight mass spectrometer (QqTOF-MS) via a CaptiveSpray ion source (Bruker Daltonics). The details of the chromatographic separation method are summarized in Supplementary Information S1 (Table S1-11). The UHPLC-QqTOF-MS/MS analysis relied on data-dependent acquisition experiments performed in the positive ion mode, comprising a survey TOF-MS scans and dependent MS/MS scans for the most abundant signals in the following 3 s (at certain t R ) with charge states ranging from 2 to 5. The mass spectrometer settings are summarized in Table S1-12.

Data Post-Processing and Analysis
The acquired raw LC-MS data were processed with PEAKS Studio software (v. 10.6, Bruker Daltonics) (for the search settings, see Table S1-13). Thereby, the identification of peptides and annotation of proteins relied on a search against amino acid sequences of Cucumis sativus cv. Chinese Long v2 [75] (Cucurbit Genomics Database v1) accomplished with the SEQUEST algorithm [99,100]. Proteomics data were post-processed in RStudio (posit.co, accessed on 10 March 2023), which is an R programming language environment [101]. Quantitative protein analysis was accomplished on log2-transformed normalized data using the limma package with the minimal abundance alteration cut-off of 1.5 logarithmic fold change [102]. Additionally, Mercator4 (v. 2.0, plabipd.de/portal/mercator4, accessed on 10 March 2023) was used for protein annotations [103,104]. The subcellular localization of proteins was defined using WoLF PSORT [105]. Enrichment analysis was carried out using KOBAS-i [106].

Statistical Analysis
All boxplots were prepared using R programming language [101] and RStudio software (posit.co, accessed on 10 March 2023). The default code for the boxplot and stripchart functions from the base R package were used. Statistical analyses of qPCR, CsRALF34 treatment, and I LRI data were performed using Wilcoxon's signed-rank test from the base R package. For stress response markers (e.g., H 2 O 2 , TBARS, 13(RS)-HPOD, Asc), the Student's t-test for independent samples combined with a normality test was applied. t-tests with Benjamini-Hochberg false discovery rate correction were also performed for the identification of metabolites and proteins up-or downregulated during CsRALF34 overexpression. In all statistical tests, differences with p-values < 0.05, <0.01, or <0.001 were considered statistically significant.

Conclusions
The results of our study on the function of the cucumber ortholog of Arabidopsis RALF34 showed that in contrast to Arabidopsis, CsRALF34 did not affect lateral root initiation at all but did affect root growth. At the proteome level, the overexpression of CsRALF34 led to increased levels of protein phosphatase 2C (PP2C) and calcium-dependent protein kinase 12 (CDPK12), i.e., regulators involved in phytohormone signaling, developmental processes, and the response to several stress factors. These effects should lead to the suppression of CDKA activity and, thus, to a block of G2/mitosis transition. Reduced levels of proliferating cell nuclear antigen and helicases are consistent with reduced mitotic activity in the basal part of the root meristem apical and might also explain the reduction in root growth. Reduced levels of NADH-dehydrogenase and succinyl CoA synthetase suggested suppression of the mitochondrial electron transport, which would be consistent with the reduced ROS levels in CsRALF34 overexpressors. However, it is difficult to determine which one has the primary effect.
Altogether, these results show that a better understanding of RALF34-receptor interactions in non-model plants should be key to understanding the physiological role of these small signaling peptides.